Nonreciprocal switching thresholds in coupled nonlinear microcavities 
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A novel concept for the design of nonlinear optical diodes is proposed which uses the mult ist ability of coupled 
nonlinear microcavities and the dependence of switching thresholds on the direction of incidence. A typical 
example of such diode can be created by combining two mirror symmetric microcavities where modes of the 
opposite parity dominate. It is shown that a strong nonreciprocal behavior can be achieved together with a 
negligible insertion loss. To describe the dynamical properties of such systems, a model based on the coupled 
mode theory is developed, and a possible implementation in the form of multilayered structures is considered. 

OCIS codes: 230.3240, 230.1150, 230.4320, 140.3945, 190.1450, 230.4170. 



The development of integrated photonic circuits requires 
ah-optical elements for high-speed processing of light sig- 
nals. The optical diode (OD) is one of such indispens- 
able elements. Similar to electronic diodes, it allows the 
flow of light only in one direction and can reduce prob- 
lems caused by unwanted reflections or interference ef- 
fects. The successful design of an OD relies on the break- 
ing of time-reversal symmetry [1]. The Faraday effect in 
magneto-optical media is mostly used for this purpose. 
There has been a considerable progress during recent 
years in making such devices suitable for integration into 
waveguides on a chip [2-4]. However, the need to ap- 
ply a strong external magnetic field still remains a ma- 
jor limiting factor for their development. Alternative ap- 
proaches to fabricate ODs are based on chiral media such 
as cholesteric liquid crystals [5] or metamaterials [6,7]. 
Although their figures of merit are significantly lower, 
they can be produced with very small sizes and do not 
require any external fields. 

Nonlinear effects attracted a lot of attention as a 
means of achieving unidirectional transmission. They 
can be subdivided into two large groups. The first one 
uses the nonreciprocal conversion of waveguide modes 
due to asymmetrically placed defects in directional cou- 
plers [8] and quasi-phase-matched gratings [9,10]. As an 
option, the asymmetry can be created by launching addi- 
tional signals that modulate the refractive index dynam- 
ically, leading to indirect photonic transitions [11]. The 
second group uses the folding of resonances in microcav- 
ities (MCs) possessing the Kerr nonlinearity and can be 
implemented in many different configurations [12]: from 
photonic crystals with embedded defects of various form 
and dimensionality [13-15] to quasiperiodic [16, 17] and 
aperiodic [18,19] structures. 

The resonant properties of the MCs allow ODs of 
smaller sizes to be constructed while keeping a large con- 
trast ratio between the transmission in the forward and 
backward directions. However, a common problem in the 
designs based on the usage of MCs is that the high con- 
trast ratio is often achieved at the expense of lowering 



the maximal transmission in the forward direction. In 
this Letter, we show how to solve this problem and to 
create nonlinear ODs with a negligible insertion loss. 

It is known that the mirror symmetric MCs demon- 
strate resonances of perfect transmission. They can be 
used as a building block to construct more complex 
structures, and if these MCs are designed to have a 
common resonance of perfect transmission, the resulting 
structure will also show it in the linear spectrum [20] . To 
create asymmetry, it is sufficient to use two such MCs 
and to modify their sensitivity to the intensity of incident 
waves. Without any loss of generality, we can consider a 
multilayered geometry shown in Fig. 1(a). In this case, 
the nonlinear response of each MC can be changed ei- 
ther directly by varying the thickness of the nonlinear 
defect region or by modifying the localization strength 
via the thickness of the Bragg mirrors. These two de- 
grees of freedom can be rigorously taken into account in 
the framework of coupled- mode theory [21,22], and we 
will investigate how to use them efficiently to preserve 
the perfect transmission in the nonlinear case. 
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Fig. 1. (Color online) (a) A schematic representation of the 
proposed OD in the form of multilayered structure. It consists 
of two MCs filled with nonlinear material, (b) A sketch of the 
model used to describe the transient dynamics. 
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The coupled- mode equations for the mirror symmetric 
MC can be written as 
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where A{t) is the amphtude of the mode Ep^{x) with the 
resonant frequency uoq and damping constant 7a, u±{s±) 
are the amphtudes of the forward and backward prop- 
agating plane waves on the left (right) side of the MC 
[Fig. 1(b)]. The Kerr nonlinearity enters via the charac- 
teristic intensity /a and is responsible for the nonlinear 
shift of resonant frequencies. If the mode E^ix) is even 
(odd), the lower (upper) signs should be used in these 
equations. 

It is convenient to go into the frequency domain and 
to derive the transfer matrix Ma which relates the plane 
waves on the opposite sides of the MC as [ix^, = 
Ma[s+, S-Y' where 
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and I is the identity matrix. Regardless of the non- 
linear term, the form of the transfer matrix (3) satis- 
fies such fundamental physical properties as the time- 
reversal symmetry (M22 ^11, M12 ^21)^ con- 
servation of energy flow (detM = 1) and in addition the 
spatial mirror symmetry (M12 = — M21). 

The transmission spectrum for waves incident from 
the left can be found as Ta = = |(Ma) 
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It follows from the formula (4) that the condition of per- 
fect transmission can be written as /in/^A = — Aco'/7a, 
where Acj = uo — uoq is the frequency detuning. This 
condition can be generalized to the case of two or more 
MCs which have the same resonant frequency but dif- 
fer in their characteristic intensities and quality factors 
Q = uJo/{2j) 

Qa/a = Qb/b. (5) 

It is useful to introduce the parameter / > 1 as 
a measure of asymmetry between these MCs / = 
Qa/Qb = Ib/Ia' Since larger quality factors generally 
imply smaller characteristic intensities, the condition (5) 
is not difficult to fulfill. 

The hysteresis of transmission for single MCs 'A' and 
'B' can be computed by using Eq. (4) and is shown in 
Fig. 2(a). For any value of the parameter / and the fre- 
quency detuning Acj, the resonances of perfect trans- 
mission are equally shifted by the Kerr nonlinearity and 
coincide in both MCs. However, the switching thresholds 



can be different, and a specific range of frequency detun- 
ings exists when only one of the MCs shows the bistable 
behavior. 

The hysteresis of transmission for coupled MCs 'A' 
and 'B' can be found similar to the linear transfer matrix 
method. After choosing some output intensity, the non- 
linear transfer matrices Ma and Mb defined according 
to Eq. (3) can be computed in a reversed order, giving 
the total transfer matrix of the structure and the input 
intensity as a result. It turns out however that the ma- 
trices do not commute in the nonlinear case, and their 
multiplication leads to a different answer depending on 
the direction of incidence [Fig. 2(b)]. The effect is maxi- 
mized when the MCs are of the opposite parity and can 
be explained in the following way. After reaching the sec- 
ond MC in the system, a part of the signal is reflected 
and increases the field in the first MC due to construc- 
tive interference. When 'A' is used as the first MC, the 
small increase causes the switching inside it, and the 
transmission through both MCs is greatly enhanced. On 
the contrary, when 'B' is used as the first MC, the small 
increase draws it away from the resonance so that the 
switching to a higher transmission state is strongly sup- 
pressed. To confirm these considerations, we performed 
a time domain simulation [Fig. 3]. 

In comparison to other approaches based on the usage 
of MCs, the proposed design of OD offers several impor- 
tant advantages. First of all, it shows not only a low in- 
sertion loss, but also it is stable for small variations of the 
input intensity and frequency detuning. This is mostly 
because the resonance of perfect transmission changes 
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Fig. 2. (Color online) (a) The hysteresis of transmission for 
single MCs 'A' and 'B'. Their quality factors Qb = Qa// 
and characteristic intensities Ib = //a are related by the 
parameter / = 2.1. (b) The hysteresis of transmission for 
coupled MCs 'A' and 'B' when signals are propagating in the 
'AB' and 'BA' directions. 
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its shape and flattens around the optimal set of parame- 
ters. Secondly, the switching thresholds for the opposite 
propagation directions are separated considerably. This 
prevents the excitation of a higher transmission branch 
on the hysteresis curve by waves moving in the blocked 
direction. At the same time, the waves moving in the 
allowed direction can cause an immediate switching to 
a perfect transmission even if they have relatively small 
intensities. 

We found that the above properties are more pro- 
nounced when the parameter / is in the following range 
2 < / < 3 and the typical contrast ratio between the 
transmission in the opposite directions varies from 10 
to 20. Although this number is not very large, the low 
insertion loss makes a series connection of such diodes 
to be very favorable so that arbitrary high transmission 
contrasts can be achieved in principle. 

The OD can be implemented in the form of a multilay- 
ered structure [cf. Fig. 1(a)]. Using layers of the quarter- 
wave optical thickness, the corresponding symbolic for- 
mula can be written as 

(HLf^D2^^(LH)^^ (LH)^^D2^^(HLf^, (6) 

where 'H' ('L') denotes layers with higher (lower) linear 
refractive index, 'D' is a defect layer made from a Kerr- 
nonlinear material, Pa{b) and mA(B) are integers. The 
wavelength used for the quarter wave condition deter- 
mines the resonance in the vicinity of which the nonre- 
ciprocal behavior can be observed. For the optimal per- 
formance, the linear refractive indices should satisfy the 
following relation (nn/^L)^^"^^ = ^d- In this case, the 
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Fig. 3. (Color online) (a) The transient dynamics of the OD 
for a continuous wave signal with the carrier frequency cj = 
LUo — 2.I7A and a staircase envelope. The intensities of the in- 
put and transmitted signals for the opposite directions of inci- 
dence are shown as a function of time, (b) The amplitudes of 
the MCs during the switching to a high transmission state, (c) 
When the input intensity is increased for the reversed direc- 
tion of incidence, the amplitudes of the MCs move away from 
the resonance condition \ A{t)\^/lA = \B{t)\'^/lA = -Auj/'yA- 



characteristic intensities will be inversely proportional to 
the width of defect regions, and the parameter / can be 
found as / = mA/^B- 

In conclusion, we proposed how to create an efficient 
OD by combining two nonlinear MCs of opposite par- 
ity. Our design shows a strong nonreciprocal behavior 
together with a negligible insertion loss and can be fab- 
ricated in a variety of configurations such as multilayered 
structures or photonic crystals with embedded defects. 

This work was supported by the German Max Planck 
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